Quantum effects on curve crossing in a Bose-Einstein condensate 
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Formation of atomic pairs by the dissociation of a molec- 
ular condensate or by inelastic collisions in an atomic con- 
densate due to a time-dependent curve crossing process is 
studied beyond the mean-field approximation. The number 
of atoms formed by the spontaneous process is described by a 
Landau-Zener formula multiplied by an exponential amplifi- 
cation factor due to quantum many-body effects. The atomic 
pairs are formed in an entangled (squeezed) state. The rate 
of stimulated processes depends on the relative phase of the 
two fields. 
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Introduction. — The concept of curve crossing was in- 
troduced in the pioneering works of Landau Q and Zener 
concerning transitions between two linear potentials. 
Curve-crossing theory has many applications in various 
fields of physics (see, e.g., Ref. ||^ and references therein), 
usually describing transitions between states whose en- 
ergies cross as a function of the time (a non-stationary 
crossing) or some coordinate (a stationary crossing). 

Recent developments in the physics of Bose-Einstein 
condensates (BEC) (see Refs. and references 

therein) demand taking account of many-body effects in a 
curve crossing. These effects, which cannot be described 
by the usual mean-field method, are treated here by a 
second-quantized formalism. Hereafter we call a process 
involving such effects a "many-body crossing" , as distinct 
from a "single-body crossing" in which a colliding pair of 
atoms is treated as a single body (a quasimolecule) . 

A major many-body effect on curve-crossing in BEC 
is Bose enhancement. In the photoassociation of atomic 
BEC this effect leads to Bose-stiniulated Raman adia- 
batic passage (see Ref. Q and references therein). An 
expression for the transition probability as a function of 
the rate of energy variation was obtained in [Q,^ by treat- 
ing the process as a curve crossing event in a nonlinear 
system. The nonlinearity may lead to a deviation of the 
transition probability from that predicted by the well- 
known Landau-Zener (LZ) formula. 

Although the association of atoms into molecules has 
been described within a mean-field approach (using cou- 
pled Gross-Pitaevskii equations) , the process of dissocia- 
tion of molecules into the atomic vacuum 

A2^A + A (1) 

cannot be described by this approach, since the coupling 



is proportional to the atomic field, the initial value of 
which is zero. Therefore, at least initially, the dissocia- 
tion must be a spontaneous process, which is neglected 
in the mean-field approach. 

The time-independent dissociation of a molecular BEC 
was considered in Ref. by using a numerical solution 
of the full quantum equations. It is shown there that 
the initial stage of the dissociation can be approximately 
described with good accuracy by analytical expressions 
obtained by treating the molecular field as a mean field, 
while the atomic field is treated as a second-quantized 
field. This method is similar to the parametric approx- 
imation used in quantum optics (see Ref. ||l^). This 
approximation is valid as long as the population of the 
molecular state is large and the effect of its depletion can 
be neglected. 

In the present paper we generalize the parametric ap- 
proximation to the case of time-dependent condensate 
energies. The energies can be varied by applying an ex- 
ternal time-dependent magnetic field or by changing the 
frequency of laser fields coupling the states. In general, 
the dissociation of a molecular condensate requires tak- 
ing into account the coupling to multiple atomic modes, 
in a manner analogous to methods of multistate curve 
crossing theory |^,|n|. Here we concentrate on the es- 
sential physics by considering only a single atomic mode, 
showing how quantum effects lead to a substantial initial 
enhancement of the dissociation and to the formation 
of atoms in squeezed states that are entangled (see also 

Similar quantum effects can also play an important 
role in the association of atomic condensates, where the 
molecular state is coupled to excited atomic trap states 
(which are initially not populated). Transitions to these 
states can substantially affect the condensate loss rate 
|p|,p^ and therefore should also be treated by a multi- 
mode model using a quantum approach similar to the 
one described below. 

The effects of Bose enhancement and the formation of 
entangled atoms should also appear in an atomic BEC, 
if inelastic collisions of the type 



A{0) + A{0)^A + A, 



(2) 



where ^(0) and A denote different internal states of the 
atoms, change the states of both colliding atoms. Such 
a process might be realized by applying an electromag- 
netic field, or two fields (as in a Raman process), in res- 
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onance with the two-atom transition, provided that the 
concurrent one-atom transition {A{0) A) is off res- 
onance with the apphed fields. The time variation of 
the energy can again be due to the Zeeman effect or to 
time-dependent electromagnetic field frequencies. An ex- 
perimental observation of Bose enhancement and entan- 
glement may be simpler in an atomic BEC than in the 
molecular one. 

General crossing. — Consider a system of coupled 
atomic A and molecular A2 fields, associated with the 
annihilation operators ip and V'm, respectively. In the 
case of one atomic and one molecular mode the system 
can be described by the Hamiltonian 



(3) 



where e (t) is a time-dependent energy of the atomic 
state, while the energy of the molecular state is chosen as 
the zero energy. The parameter (7am describes the cou- 
pling of the atomic and molecular fields due to hyperfine 
interaction or laser fields. 

The parametric approximation used in Ref. |^ for the 
case of constant e can be generalized to the case of a 
time-dependent e, by retaining the assumption of a large 
molecular population. We can then replace the opera- 
tor ipm by its expectation value {ipm), and obtain the 
equation of motion for the atomic field operator -0: 



(using units in which h ~ \), where 

9 = 5am (■0m)*- 



(4) 



(5) 



Let us also consider two atomic states, A{Q) and A, 
associated with the field operators V'o and '0, respectively. 
Taking into account one mode for each of the states, one 
can describe the system by the following Hamiltonian 



H = e{t) V-V + gaaV'JV'JV'V' + glJHHo^Po 



(6) 



where e [t) is the time-dependent energy of state A, while 
the energy of state A{Q) is chosen as the zero energy. 
Whenever both wavevectors fco and fc, of the entrance 
and exit reaction channels (^, are in the near-threshold 
region, the inelastic scattering amplitude should have 
the form T — \/kokb, where b is an energy-independent 
length (see Ref. |16 ). The parameter t/aa can then be ex- 
pressed as (7aa — 47r6/m. The parametric approximation 
can be applied here if the state A{0) has a large pop- 
ulation, replacing the product V'oV'o by its expectation 
value {ipoipo). If the atoms are initially in a condensate 
(coherent) state, then (ipoi'o) = (V^o)^- Denoting 

= 5aa(V'O'0o)*, (7) 

one obtains again the same equation of motion (^ for 
the field operator ip. 



The set of two first-order equations [Eq. (y) and its her- 
mitian conjugate] can be transformed to a single second- 
order equation 



(8) 



Although it is an operator equation, it is linear and its 
solution can be expressed as 



{t) = tP, (t) ^ {to) + tps (t) (to) 



(9) 



where tpcs (t) are two independent c- number solutions of 
Eq. (|). Equations (1) and (|) at i = to give the initial 
conditions for t/jcs (^o) and their derivatives: 

i^cito)^!, 0.(to) = O, (10) 
^c{to) = -ie{to), ^s{to)^-2ig*. (11) 

The solutions ?/'c {t) and ips [t] can be expressed in 
terms of a pair of conventionally selected independent 
solutions of Eq. (||) ipi (t) and ip2 (t), treated as "stan- 
dard solutions" , in the form 



V'c it) = 



(e (io)^2 (to)-«^2 {to))yJi it) 



(t) 



W^Wl,V'2} 

-(e (to) 01 (to)-i0i (io))02 it) 
2ig* 



(12) 



T^{01,^2} 



[^2 (io) 01 (t)- 01(^0)^-2(0], (13) 



where M^{'0i,'02} = ■0i0'2 — '02'0i is the Wronskian of 
the standard solutions. In the special case of a time- 
independent e 1^ , 0s {t) and Tpc (t) can be expressed in 
terms of sine and cosine functions (circular or hyperbolic, 
depending on the ratio of g to e). 

The use of Eq. (||) allows us to express the expecta- 
tion value of any function of 0^ and ip in terms of the 
expectation value at a particular time t — to- Thus, the 
number of atoms is a sum of two terms (0^ (t) -0 {t)) — 
nsp it) + ngt it) , where the first term 



Usp (t) = 1-08 it) 



(14) 



corresponds to spontaneous transitions into the vacuum 
of -0, and the second term 



nst it) = (10c it) P + 10, it) p) (^t (tg) ^ (ip)) 
+2Re{ij:it)i^c it) (0 ito)4^ito))] 



(15) 



depicts stimulated transitions taking place when the final 
atomic state is initially populated. The first and second 
terms on the right-hand side of Eq. (^5|) correspond to 
a noncoherent and a coherent initial occupation, respec- 
tively. Consider also the expectation value of a general- 
ized coordinate operator 



Q(t) =0(t)e'''+V;t(t)e"'« 



(16) 



2 



frequently used in the theory of squeezed states (see, 
e. g., 0). We present here only an expression for 
the spontaneous process, in which (to) 'tl' (to)) — 
(V- (to)V^(to)) = 0: 



(17) 



Linear crossing. — Single-body curve-crossing prob- 
lems are commonly solved by using the LZ formula, which 
is an exact solution of the linear non-stationary prob- 
lem. Let us consider a generalization to our case of a 
many-body crossing, choosing e (t) as the linear function 
e (t) — I3t. The second-order equation (||) attains in this 
case the form of a parabolic cylinder equation (see Ref. 

PI) 



i:+{l3^f-A\g\^ + il3)^j = Q. 



(18) 



This equation differs from the corresponding c-numbcr 
equation for a single-body crossing (see Ref. |]ll|) by the 
sign before A\g\'^ due to the effects of secondary quan- 
tization. An appropriate choice of a pair of standard 
solutions is then 



(19) 



^2 it) = -VW2e- 
9 



t/4. 



where U (a, x) is the parabolic cylinder function, defined 
by Eq. (19.3.1) of Ref. (l|. Here r = and 



A = 2|5|V/3. 



(20) 



The parameter A is nothing else but the LZ exponent, 
since the coupling of states, according to Eq. (|^), is 2g 
and the slope of the two-atom energy is 2/3. 

The solutions ( p^ ) have the Wronskian ^^{-01, -02} = 
i{P/g) exp(— 7rA/2) and in addition obey the conditions 



- = 2g>3% (« = 1,2), 
which allow us to express ipc (t) as 

,ri{to)^l (i) - V-S (io) 02 (t) 



0e it) = -2ig* 



W^{01,02} 



(21) 



(22) 



Using a common procedure of curve crossing theories, 
we can start by considering transitions during a finite 
time interval [t^, t], with < < i, assuming sufficiently 
large values of |io| and t (obeying the validity criteria pre- 
sented below). Asymptotic expansions of the parabolic 
cylinder functions p8[ give us at ^ 



01 [h) 



1 



■ exp 



^lx + ^\+^S{\r,\) 



V'2 (ii 



„)^iy|exp(-J 



(23) 



A-z- 



iS 



(|ro|)) 



where tq = %/2^io and S (r) = t^/4 — A Inr. The asymp- 
totic expansions at t — > oo, obtained by using the relation 
(19.4.6) in |jl||, are given by 

V-i {t) - - exp ( — A - I— iS (r) 

+ -P7 \//3 sinh (ttA) exp f ^A ~ i— ~ iS (r) — i argF («A)) 
La \4 2 / 



(24) 



V^2 (i) 



*^-*'S'(t) 



-f V^2 sinh (ttA) exp f-A - «- + iS {t) + iargT (iX) 
\g\T V4 2 

Spontaneous process. — Substituting the asymptotic 
expansions, Eqs. ( p3|) and (|2^), at large |to| and t into 
Eqs. ([l3| ) and (|lj), one can obtain the number of atoms 
formed by the spontaneous process 



n,p (t) ^ (e^-^ - 1) + 2e^yA (e2-^ - 1) 



cosx(t) cosxdrol) 



To 



, (25) 



where x (t) = /2 - 2Alnr -f argP (lA). 

Equation (|2^) contains two terms of an asymptotic ex- 
pansion in inverse powers of r and tq. The asymptotic 
approach is applicable if the second term in Eq. ( p5|) is 
small compared to the first one, or 



t,\to\ » \g/p\, 



(26) 



in agreement with estimates of the transition region for 
two-state single-body crossings (see Ref. |11 ). 

Let us recall that all the analysis above is based on the 
assumption that the depletion of the initial (atomic or 
molecular) condensate state is negligibly small, applica- 
ble when Usp (t) <C N, or 



A< — InTV, 



(27) 



where N is the initial number of atoms. Since ^ 1 (in 
current experiments an atomic BEC contains 10^ to 10^ 
atoms) the present approach is applicable even for large 
values of A corresponding to adiabatic transitions. 
The leading term of Eq. (E3) can be written as 



(28) 



The first factor here is the familiar LZ probability, while 
the second one describes an amplification of the sponta- 
neous process. This amplification is reminiscent of las- 
ing, as the exponent can be interpreted as a product of a 
characteristic crossing time (g/P) and an "amplification 
coefficient" g. 

The LZ exponent A is itself enhanced, being propor- 
tional to A^ fsee Eqs. (|)and (|2^)]. This is also true for 
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the association process, but unlike that process, the Bose 
enhancement introduced here by the amphfication factor 
is unique to the dissociation process. As an example, 
consider the Na Feshbach resonances studied in Ref . . 
The parameter A there can be of the order of unity, de- 
pending on the ramp speed, in which case one obtains an 
amplification factor of more than a 100. 

A distinction between many-body and single-body 
crossings appears also in the interpretation of the cross- 
ing as scattering on a parabolic potential barrier. Indeed, 
on replacing the time by a coordinate, Eq. (|^) is trans- 
formed into a corresponding one-dimensional stationary 
Schrodinger equation. In the single-body case the energy 
lies above the barrier (see Ref. |l^]). However, in the 
many-body case, as one can see from Eq. (H), the energy 
lies below the the barrier, and the process involves bar- 
rier penetration. This leads to the replacement of the 
negative exponents by positive ones, and yields the first 
term of Eq. ( |25|) in place of the LZ formula. 

Using Eqs!^l|), (||), (|2|), and (||), we can write Eq. 
©as 



y/e^""^ - 1 



+ e^p{TrX-ixiT) -iargg + 2i0-iTT/2) . (29) 
Substituting one of the two orthogonal phase angles 9 = 



Eq. (^] or the atomic field (ipo) [see Eq. (Q)], the stimu- 
lated dissociation is indeed dependent on the phase dif- 
ference between the two fields (a and (V'm) or {ipo))- This 
effect is reminiscent of the amplification or absorption of 
coherent light in laser-active media, and its dependence 
on the phase of the light. 

Conclusions. — The parametric approximation al- 
lows us to obtain an analytical expression (^^ for the 
number of atoms spontaneously formed by the crossing 
of linearly time-dependent energies of the atomic and 
the molecular states or two atomic states. Many-body 
effects lead to an enhancement factor in the LZ expo- 
nent, as well as to the multiplication of the LZ proba- 
bility by the amplification factor in Eq. (|2|) describing 
a positive-exponential enhancement of the spontaneous 
process. The atoms are formed in a squeezed state. The 
number of atoms formed by the stimulated process de- 
pends on the relative initial phase of the two fields. Sim- 
ilar quantum effects should also apply in other cases in 
which particles are created in pairs. An experimental 
realization of the processes we discuss may allow produc- 
tion of an entangled atomic gas. 

The authors are grateful to P. D. Drummond and K. 
V. Kheruntsyan for helpful discussions. This work was 
partially supported by the NSF through a grant for the 
ITAMP at Harvard University and the SAO. 



d± = h{x{r) + argg) ± f in Eq. (|6|) we obtain 



{Q±)sp 



Ve 



2tt\ 



A>1 



2ttX 



(30) 



Therefore, the spontaneous dissociation forms atoms in 
a squeezed state (similar to a squeezed state of light 1 17 ) 
in which has an increased uncertainty and has 
a reduced uncertainty. A similar effect concerning the 
time- independent case has been described in Ref. Q . 

Stimulated process. — Consider now the transition 
stimulated by a finite initial occupation of the final 
atomic state. Starting from a coherent state \a) of the 
operator ■(/', such that t/j (tg) |a) = a\a), and using Eqs. 
(|l|), (H), and (H), one obtains 



v/e^'r^ - 1 



-l-exp (ttA -f ix (to) -f iarg5 -I- 2i arga -f i7r/2) . (31) 

It follows from this expression that the final atomic pop- 
ulation is dependent on the phase of a, lying in a range 
bounded by the maximal and minimal values 



n± 



Ve 



2itX 



lie 



2ttX 
,-2Tr\ 



(32) 



at arga — —^{xihoD + argg) =F |-, respectively. Since 
argg contains the phase of the molecular field (V'm) [see 
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